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Abstract 

We first study Zn walls/interfaces in a deconfined phase of Witten's DA- 
brane background of pure SU(N) Yang-Mills theory, motivated by a recent 
work in the case of N = 4 SYM. Similarly to it, we propose that for a large 
wall charge k ~ N, it is described by k D2-branes blown up into a NS5- 
brane wrapping S 3 inside S" 4 via Myers effect, and we calculate the tension by 
suitable U-duality. We find a precise Casimir scaling for the tension formula. 
We then study the fate of -Z^-vacua in a presence of fundamental flavors in 
quenched approximation via gauge/gravity correspondence. In the case of 
D3/D7 system where one can vary the mass m q of flavors, we show that there 
is a phase transition at T c ~ m q , below which the Z^- vacua survive while they 
are lifted above the critical temperature. We analytically calculate the energy 
lift of fc'th vacua in the massless case, both in the D3/D7 system and in the 
Sakai-Sugimoto model. 
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1 Introduction 



Gauge/gravity correspondence has often provided us easy tools for addressing difficult 
problems of strongly interacting gauge theories. It is based on the observation that large 
N, large t'Hooft coupling limit of certain gauge theories have alternative weakly-coupled 
effective descriptions in terms of gravity or string theory One important aspect of 
these dual theories is that the spacetime they are living is one-dimension higher than 
the original spacetime of gauge theories, and the additional dimension, often referred to 
as holographic coordinate, plays a physical role of energy scale of a given process. This 
mapping of space-energy is practically achieved by the presence of a warping factor in the 
metric that depends on the additional dimension, which provides a geometrical potential 
along the holographic dimension. 

Out of vast amount of previous works that have been done over a decade to test and 
study this idea, one interesting application has been to describing the deconfined phase 
of gauge theories in a sufficiently high temperature, whose dual picture is a black-hole 
spacetime with its horizon located at some point in the holographic coordinate [2, [3] . Due 
to the geometric barrier along holographic dimension provided by the warping factor in 
the metric, this black-hole spacetime is stable against its Hawking radiation, resulting 
in an equilibrium state called Hartle-Hawking state. Performing relatively easy analyses 
in this black-hole background has given us many valuable, typically non-trivial, results 
about physics of high temperature, deconfined gauge theories in strong coupling [3]. 

In this work, we would like to add one more piece to the mounted pile of results about 
high temperature phase of gauge theories obtained using gauge/gravity correspondence. 
Our motivation is two-fold. We first study Zn walls in the Euclidean picture of finite 
temperature deconfined phase of pure SU (N) Yang-Mills theory using the Witten's dual 
background of ZM-branes [3]. We are interested in the /c-wall with 4 ~ which is 

largely motivated by a recent work of Armoni, Kumar and Ridgway for the case of N = 4 
SYM via the AdS^ x S 5 background [5] . Similar to their result, we will show that the A;- wall 
with jj ~ 0(1) is k (Euclidean) L>2-branes blown up to a single iVS'5-brane wrapping S 3 
inside the S 4 in the background via Myers effect. We calculate its tension after performing 
suitable U-dualities of Type II SUGRA to make the problem more tractable. Our result 
of fc-wall tension indicates certain limitations of finite temperature black-hole geometry 
of Witten's background that has been used in the previous literature. 

Our second motivation is to study the fate of Zn discrete vacua, and hence Zn walls 
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connecting them, when we include fundamental flavor quarks in the gauge theory. One 
naturally expects that since Zn is no longer a symmetry of the theory in the presence 
of fundamental matters, these Zn vacua are generally lifted except the one true vac- 
uum. However, there appears an issue when the fundamental matters are much more 
massive than the temperature, because the physics with energy scale below the mass of 
fundamental matters decouples from the fundamental matters, and the effective low en- 
ergy theory would still have a center Zn symmetry. One logical possibility is that when 
■^r ^> 1, the Zn vacua and domain walls persist in the theory, while they disappear at 
-C 1, and these two phases are separated by a phase transition at ^ ~ 0(1). We will 
give an evidence to this scenario in the framework of gauge/gravity correspondence using 
-D3/-D7-brane system corresponding to N = 2 fundamental hypermultiplets in N = 4 
SYM ©. 

In the case of massless fundamental flavors, Zn vacua are lifted unambiguously. Using 
gauge/gravity correspondence, we calculate analytically this energy lift of k- vacua in 
quenched approximation, in both D3/D7 system and in the Sakai-Sugimoto model [7]. 
As a by-product, we also find an infinite tower of stable states above each lifted k- vacua, 
similar to the tower of stable states above vacuum that was found in Ref.[8]. 

Note added in revision: 

The author was informed of an issue of correct interpretation of center Zn symmetry 
and the corresponding Zn walls in Euclidean thermal SU(N) gauge theory 3- If one takes 
SU (N) I Zn as the gauge group rather than SU (N) , the center Zn symmetry is absent and 
the Zn walls are not the Minkowski domain walls, but instead some Euclidean instanton- 
like objects [11]. In the presence of fundamental flavors, the gauge group is indeed SU(N) 
and the center Zn symmetry can be at least an approximate symmetry of the theory. Our 
second subject of the paper should be viewed in this way. Alternatively, one can think of 
thermal S 1 as a spatial compactification, and consider 3-dimensional Minkowski theory 
after Wick rotation along one non-compact spatial direction. In that context, one can 
talk about Minkowski domain walls in the 3 dimensional theory. 



1 We thank Andrei Smilga for bringing this issue to our attention. 
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2 Review of center Zn in gauge/gravity correspon- 
dence 

This section is devoted to a brief review on the center Zn symmetry of the Euclidean 
description of finite temperature SU(N) gauge theory with only adjoint matters, and its 
spontaneous breaking in the deconfined phase via non-vanishing Polyakov line along the 
thermal circle 0, [TQ] . We will also summarize relevant results in Refs. [SI [12] in the context 
of gauge/gravity correspondence to lay ground for our analyses in subsequent sections. 
Readers familiar to these subjects can jump into the next section. 

2.1 Center and its spontaneous breaking in hot SU(N) gauge 



An Euclidean finite temperature gauge theory is defined in the space S 1 x i? 3 , where the 
thermal circle, or Euclidean time, has a period (5 = A. As a SU(N) gauge theory, field 
configurations that are connected by SU(N) gauge transformations must be treated as 
identical to each other, and one should mod out the space of fields by gauge transfor- 
mations. Having S 1 means that the local gauge function U(x) G SU(N) that we mod 
out should be periodic along S 1 to preserve usual periodicity of matter fields that are 
charged under SU(N). If the theory contains only adjoint matters, one might think of an 
extended notion of gauge transformation whose gauge function U (x) G SU (N) is periodic 
only up to an element of the center inside SU(N). This is because for adjoint matters 
$, as Z N freely commutes with adjoint $, the gauge transformation 



doesn't change the periodicity of $, and seems to be allowed without any problem. Note 
however that under the extended gauge transformation, the Wilson line along the thermal 
S 1 , called Polyakov line, 



transforms exactly by the center Z N element W{S 1 ) — > e J N~W(S 1 ). Therefore, if we 
accept WIS 1 ) as one of physical "gauge invariant" observables of the theory, the extended 
gauge transformations are in fact not allowed as gauge transformations that we mod out. 
Rather they should be thought of as symmetries of the theory as the action is clearly 



theory 



-> U(x)*${x)U(x) 



(2.1) 




(2.2) 
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invariant under the transformations. The extended gauge transformations are constant 
Zn times the usual periodic gauge transformations that we mod out, so that the resulting 
symmetry is a global symmetry. 

One can immediately identify WIS 1 ) as an order parameter of the spontaneous break- 
ing of this center Z^ symmetry. In fact, it can be argued to get a non- vanishing expec- 
tation value in the deconfmed phase of high temperature, while it should vanish in the 
confined phase. The rough picture is that WIS 1 ) represents a world-line of an external 
fundamental quarks sitting at a point in R 3 in Euclidean finite temperature description. 
In the confined phase, its presence would cost too much free energy due to confinement so 
that the partition function with WIS 1 ) would vanish, while one generally expects a finite 
non-vanishing result in the deconfmed phase. Therefore, Z^ is spontaneously broken in 
the deconfmed phase, and there exist iV-number of vacua parameterized by the Z^ phase 
of the order parameter W{S 1 ). One naturally thinks of domain walls connecting these 
vacua. By fc-wall, we will refer to a domain wall connecting i'th vacuum and % + fc'th 
vacuum, and k is identified with k + N by definition. 

There is a beautiful connection between a Z^ domain wall and a large spatial t'Hooft 
line |13j . Consider a large spatial t'Hooft line along a curve C which bounds a large 
surface S spanning two space directions in R 3 , say x 1 and x 2 . By t'Hooft line, we mean 
inserting a magnetic monopole whose world-line is the curve C. In the presence of a 
magnetic monopole, the Bianchi identity for the electric gauge potential A is violated by 

dF = 5f , (2.3) 

(3) 

where 5 C is the Poincare dual 3-form to C. Being familiar to the case of Dirac monopole, 
we can think of S as a Dirac sheet, a world-sheet spanned by a Dirac string. As dS = C, 
we have 

dSf ) = 5§ ) , (2.4) 

and we see that 

F = dA = 4 2) • (2.5) 

Then, consider a cylinder D made of S 1 times a finite interval in R 3 which connects two 
points P, Q from the two opposite sides of the surface S, so that this interval crosses 5* at 
a single point in R 3 . It is clear that this cylinder has an intersection number 1 with the 
surface S in our spacetime S 1 x R 3 , or equivalently 

/ 4 2) = #(£>, S) = l . (2.6) 

JD 
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Using (12.51) . the left-hand side in the above is 



F= A= A- A , (2.7) 

D JdD JS 1 atP JS 1 atQ 

which tells us that the phase of Polyakov line jumps across the surface S of Dirac sheet, 
and we can think of S as a domain wall separating regions of different Polyakov lines. 
Although the above discussion is given in terms of an Abelian gauge theory, the logic is 
essentially identical in the case of non- Abelian SU (N) gauge theory, and the conclusion 
is that the /c-wall can be thought of as the minimal surface S bounded by a large spatial 
t'Hooft line of k-th anti-symmetric representation of the magnetic group. 

2.2 Zn vacua and Zjy walls in gauge/gravity correspondence 

As one typically considers the large N limit in discussing gauge/gravity correspondence, 
it may at first seem unlikely to see discrete symmetry in the gravity dual descrip- 
tion. However, believing AdS/CFT correspondence beyond the leading large N limit, one 
should be able to access sub-leading effects by carefully taking into account quantum 
effects in the gravity side, which is not easy in general. In the case of Z^, Aharony and 
Witten successfully identified the relevant quantum effect in AdS§ x S 5 background, that 
is responsible for the discrete Zn symmetry appearing in the gravity dual description of 
the N = 4 SYM [121. 



For our purposes, we will recapitulate only the case of Poincare patch version of 
AdS§ x S 5 corresponding to the gauge theory defined on S l x R 3 , while we refer the 
readers to Ref.[12| for the discussion of N = 4 SYM theory on S 1 x S 3 . This is because 
only in the former case, we can have a spontaneous breaking of global symmetry such as 
Z^ that we are interested in. A finite temperature phase of N = 4 SYM in Euclidean 
description has a unique dual geometry given by the Euclideanized black-hole in the 
Poincare patch; 



ds 2 - — 
as E - ^ 2 



' vr 4 T 4 L 8 \ , 2 , 2 
1 — I dt E + dx { 



where the thermal circle combined with the holographic radial coordinate r makes a two- 
dimensional cigar shape D that closes off at the location of the horizon r# = nTL 2 with 
L 4 = 47rg s iV7 4 . The theory is always in the deconfined phase due to conformal symmetry. 
As is well-known, the expectation value of a Wilson line is calculated by the semi-classical 
string world-sheet which has a boundary at UV r — > oo along the curve of the Wilson line 
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[IH[T5]. O ne easily identifies such string world-sheet in the case of our Polyakov temporal 
Wilson line WIS 1 ); it spans precisely the cigar D of the thermal circle and r. Because the 
world-sheet closes off at a finite point r = r#, it would give us a finite value of the string 
action, and hence non-vanishing expectation value of WIS 1 ) after suitable holographic 
renormalization [16J, which is in accord with the fact that W(S X ) is non-vanishing in a 
deconfined phase. However, there appears a puzzle since it seems one can freely turn on 
NS 2-form B with dB = on the cigar without any change in the Type IIB SUGRA 
action, with arbitrary values of 

2^L B • < 2 ' 9 > 
which appears as a phase of our semi-classical string amplitude, and hence of W^S 1 ). As 
this B mode is a normalizable mode, we in fact have to sum over these possible phases 
in the partition function, which would make the expectation value of WIS 1 ) zero [3]. 
To prevent this, there must be a mechanism that lifts the degeneracy among continuous 
values of NS B field, and we indeed also need this to have discrete values of the phase 
of WiS 1 ). 

Ref . [T2] showed that the necessary mechanism is from quantization of the RR 2-form 
C in our background with N D3 flux on S 5 , 

' F 5 = N , (2.10) 



(2vr/ s ) 4 Js* 
such that the action contains the piece 

N- ( -j-r / B ) • . \ „ / dC = 6 ■ \ N0 / dC , (2.11) 
\2h1 2 Jd j (27r/ s ) 2 ii?3 (2-nl s ) 2 jR? ' v ; 

from the Type IIB term J Mw F 5 A B 2 A H 3 with H = dC. With the usual kinetic 

term for H = dC along R 3 obtained after integrating out S 5 x D, 



9sN f .„ N2 _ 1 



2V/*T3 Jm iHl23r = 2? Ls {Hl23)2 ' (2 ' 12) 

the above 9 was shown to play a similar role as the ^-F i term in 2D Abelian U(l) theory, 
where the vacuum energy is given by 

E e = ■ mm kez (2nk - 6) 2 , (2.13) 

with integer k labelling certain allowed quantum states of the system. This says 9 = 2ttZ 
are discrete vacua among continuous possibilities, and looking back (12.111) . one readily 
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finds discrete Zn vacua of 

1 r 1-nk 

™iL B =ir • k ~ k+N • < 2 - 14 > 

in the quantum theory, which at the same time also explains the Zn phase of W{S 1 ) at 
each vacua. 

The Z^ walls in the gravity description can most easily be identified by considering 
their connection to large spatial t'Hooft lines that we discussed above. In the simplest 
case of k = 1 wall, the t'Hooft line corresponds to a (Euclidean) Dl brane world-sheet 
whose boundary at r — ► oo is along the t'Hooft line curve. One may think of the Dl 
brane as the Dirac sheet of the t'Hooft line, and one concludes that Dl brane spanning 
along two spatial directions in R 3 is the k = 1 domain wall. Note that when the t'Hooft 
line is large, most of the Dl world-sheet will sit at the IR end r = r# to minimize its 
tension. Therefore, the Dl brane sitting at r = r# and a point in S 5 , while spanning two 
directions in R 3 is the final configuration we are seeking for the k = 1 wall in N = 4 SYM. 
In a recent work by Armoni, Kumar and Ridgway [5], they proposed that for fc-walls with 
jj ~ 0(1), the k Dl branes are blown up to a NS5 brane wrapping S* 4 C S 5 via Myers 
effect, and the result for the A;- wall tension in their picture supports this claim, that is, 
for low k the tension is simply that of the k Dl branes, and (N — A;)- wall tension is equal 
to the fe-wall tension. 



3 /c-walls in the Witten's I}4-brane background 

We come to one of our objectives in this work; calculating fc-wall tension with ~ 0(1) in 
the (approximate) gravity dual of pure SU(N) Yang-Mills theory. We first discuss briefly 
how center Zn appears in the Witten's background, which is rather closely parallel to the 
case of AdSs x S 5 . 

There are two finite temperature Euclidean geometries competing with each other 
near the confinement-deconfinement phase transition; a confined phase at low T where 
the thermal circle remains finite in all spacetime region, and a deconfined phase of Eu- 
clideanized black-hole where the thermal circle shrinks to zero at the horizon, making a 
cigar shape combined with the holographic direction [T7] . 

The confined phase geometry is 



7 



^ - 2 



-e 4 , e v = g s — , V4 = Vol(5 ) = 



V 4 Vi? 
i? 3 = ixg s Nl 3 s , /(f/) = l-(^) 3 , (3.15) 

with the thermal circle tE periodic ~ ig + /3, and 24 is compactified with the period 

The useful relation between the above parameters and the 4D gauge theory data on (tE, %i) 
is given by [7] 

R 3 = 9 -lw^ > ^ = ^ M iVM^ , ^=^y ■ (3.17) 
2M X ic 9 lnM KK l s 

As usual the Polyakov line expectation value is described by a string world-sheet wrapping 
tE and one more direction in the geometry, and since tE never shrinks to zero one easily 
finds that there is no way this string stops in the geometry without turning back to the 
UV boundary U — *■ 00, which corresponds to another anti-Polyakov line in the gauge 
theory. In short, (W(S 1 )) = with a single insertion of W(S 1 ) in the gauge theory, in 
accord to the expectation in a confined phase. The way discrete Zn symmetry appears in 
the situation is that if we insert W(S 1 ) , s by multiple N times, these multiple N number 
of string world-sheets can now sit on (Euclidean) DA branes wrapping 5 4 x because 
on this baryonic D4-brane, there is a Chern-Simons coupling 



/ F A A (2nl 2 s )A = — -i— [ F,AA = N f dt E A , (3.18) 

J DA (2lll s ) 6 JDA J 

which induces N string charges on the world- volume, so that N strings can sit on it [18] . 
Therefore, multiple N times insertion of W / (5' 1 ) , s gives us a finite expectation value and 
one concludes that there is a Zn symmetry under which W(S 1 ) has a unit charge. Note 
that the Z^ symmetry is unbroken in this phase due to (W(S 1 )) = 0. 

Our present interest concerns more about the deconfined phase at high T, for which 
the dual geometry is given by 
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- ( , , 1 , -A /R\2 I dU 2 



v 4 a \R 

R 3 = 7rg s Nl 3 s , f(U) = i-(!k) , (3.19) 
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where the period of the thermal circle, (3 = k, is related to Ut by 

n An fR 3 \^ , 

and we have re-scaled x^ such that its period is (2nl s ) for later convenience. In this 
background, the string world-sheet for a Polyakov line can stop at U = Ut where the 
thermal circle shrinks to zero making a cigar shape D with U > Ut, so that the string 
action becomes finite and one expects a non- vanishing VEV of W(S l ). The same issue 
of NS 2-form phase Id B ns arises as in the case of AdS§ x S 5 , and the resolution is 
also similar. In Type IIA, there is a term 

J^tf Lo F 4 A F 4 A 5 2 = iV ■ [7^2 J D B NS j F^e- — F 4 , 

(3.21) 

and considering quantizing F± R along (xi, X4) after integrating over S 4 x D to get the 
kinetic term 

3g a N 

2 9 7T 6 / S 5 T 3 

we have the integer parameterized vacua of 9 = 2nZ to minimize the energy of the allowed 
quantum states treating x\ as a time, 



J d Xl dx A (F1234) 2 , (3.22) 



» t = «P U J a.™) = exp ( ™ / q«) . (3.23) 



This provides us the Z^r vacua of 



' r B NS = ^ , A;~A; + iV , (3.24) 



2nl 2 s AT 

corresponding to the spontaneous breaking of Zn by (W(S 1 )) ^ 0. 

It is not difficult to identify the string theory object that plays a role of do- 
main walls, at least for fc-walls with sufficiently low fc's. Again, using the large spatial 
t'Hooft line is most convenient. In our DA background, the magnetically charged sources 
compared to the electrical degrees of freedom of fundamental strings are provided by £4 
wrapping D2-branes [T9]. This can also be seen by considering T-duality along x 4 upon 
which N Z)4-branes become D3-branes and D2 becomes Dl, the usual magnetic degrees 
of freedom. Consequently, a large spatial t'Hooft line will be the boundary of the world- 
volume of a a^-wrapped (Euclidean) D2-brane at U — > 00, and most of the D2-brane 
world-volume will sit at the IR end U = Ut to minimize its tension. We can think of 
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two spatial directions in R 3 that D2-brane is spanning as the Dirac sheet corresponding 
to the t'Hooft line, and hence the k = 1 domain wall. Therefore, the D2-brane sitting at 
U = Ut and a point in S 4 while spanning x 4 and two directions in R 3 is the desired k — 1 
domain wall in the gravity side. For small number of k, one naturally expects a collection 
of k -D2-branes to be the corresponding object for the fc-wall. Its tension is calculated 
straightforwardly by the DBI action of D2-branes at U = Ut, 



U=U' 



This result is qualitatively different from both the weak-coupling calculations of pure 
SU(N) Yang-Mills theory [201 El]> and the recent weak/strong coupling calculations of 
N = 4 SYM by Armoni, Kumar and Ridgway |5J, where the results have a common form 
of 

T 2 

T k = (const) ■ kN , (3.26) 

for small <C 1. 

One can understand the origin of the discrepancy as follows. The confinement- 
deconfmement phase transition in the Witten's background happens at a temperature 
T c of order Mkk P2j, and the above deconfinement background becomes relevant when 
T is much larger than Mkk- However, Mkk also serves as a compactification scale of X4 
below which we have a 4D gauge theory, while above which massive KK modes enter the 
dynamics and the theory becomes effectively 5-dimensional. The number of KK modes 
that would enter at the temperature T is roughly given by m kk 1 and assuming each 
new degrees of freedom contributes to the tension of the domain wall, one can understand 
= ' Mkk behavior of our result for the tension formula. However, missing \J gy M N 
factor remains still puzzling. One possibility is that the mass of KK modes might become 
7== at strong coupling, and the number of effective degrees of freedom at T might 

V9ym N 

be m^kk \l Oym-N instead of M T . Our result may be considered as pointing out this phe- 
nomenon. The above discussion indicates that the deconfinement geometry (I3.19P that 
has been used in the previous literature has some limitation to be used as a gravity dual 
background of a 4D gauge theory in its deconfined phase. 

We proceed by studying /c-walls with k being large and comparable to N, motivated 
by a recent work of Armoni, Kumar and Ridgway in AdS§ x S 5 [5J, and our method in 
this regard will be similar to theirs. The naive picture of fc-wall as a simple collection 
of k .02-branes should break down when k ~ iV because k is Zjv-valued and moreover 
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k is related to (N — k) by charge conjugation and the tension must be invariant under 
this. Similar to the claim in Ref . [5] , we propose that k D2-branes sitting at a point in S 4 
blow up into a single IIA NS 5-bran.e wrapping S 3 C S 4 due to the background _F 4 flux on 
S 4 , and we confirm this picture by computing the fc-wall tension and check the necessary 
properties. 

To analyze more easily the Type IIA A^S"5-brane dynamics, with k /}2-brane charges 
on its world- volume along 24 and two spatial directions in R 3 , say x\ and 22, we perform 
U-dualities of Type II theories in the following way. Note first that in the deconfined 
phase geometry (13.191) . 24 size remains finite in all region of the spacetime, contrary to 
the confined phase geometry where 24 shrinks zero at U — Ukk- Hence, we are eligible 
to take T-duality along 24 in the deconfined phase. We stress that this is not allowed in 
the original Witten's background with shrinking 24, and one should take this operation 
applicable only to the present deconfined phase as a mere tool for calculating some physical 
quantities in an easier way. After T-duality, our A^S"5-brane becomes Type IIB iVS'S-brane 
because it wraps the 24 direction, and the k D2-brane charges wrapping 24 will transform 
to k Dl-brane charges spanning now only Xi and 22. They are homogeneously distributed 
in 24 and S 3 C S 4 of our A^S"5-brane world-volume. We next perform Type IIB S-duality 
in the system, such that we finally have a (Euclidean) D5-brane wrapping S 3 C S 4 and 
s 4,i,2) with k Fl-string charges on its world-volume along 2^2, homogeneously distributed 
in 24 and S 3 C S 4 . This can be studied by DBI plus Chern-Simons action. 

The T-dualized background of (13.191) along 24 looks as 



3 3 

d~s 2 = (|) , (/(^ + ^) + (|)^^0 2 ^+(§)*(^ + ^dn 

(2vr/ s ) 3 iV 



d2 4 Ae 4 , e* = e+—= = g.(M KK l a ) , (3.27) 



where the parameter R and f(U) are same as before, and more importantly 24 has the 
same period (2nl s ) as 24. One should remember that the relevant rules of T-duality in 
Ref. 1 22 1, that is, 044 = — and = e^—^= must be applied in the coordinate with the fixed 



n f 44 vy 44 

period (2ttI s ]q. The resulting 5-form flux simply describes iV D3-branes homogeneously 

distributed along 24, and the dilation is constant as it should be in a D3 background. 

After a further S-duality, upon which ds' 2 = e~^ds 2 and e^' = e~^, the final geometry is 



2 Moreover, C^ c = 4C*| hcrc and C$ C1C = |C^ horo . 
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given by 



g s {M KK l s ) \RJ \ > g s \U J 



g s (M KK l s ) \UJ \f(U) 



dU 




+ 

9s[M 

F 5 = ^£ 4 A £4 , e? = ; , (3.28) 

V* g s {M KK l s ) 

and we consider a D5-brane wrapping S 3 C S 4 and £4,1,2 with k Fl charges along x\^- 
The (Euclidean) D5-brane action is (we omit primes in the above geometry) 

= ^J d 6 £ e'ty/det (g* + 2nP s F) - z> 5 J Cf R A 2tt1 2 s F , (3.29) 

with n v = (2tt)-p/ s -(p +1 ) and we turn on the world-volume gauge flux F = dA only along 
£1,2 to represent k Fl charges. To find an expression of C RR usable for our purpose that 
satisfies dC RR = F 5 with F 5 given above, we introduce a polar angle < 9 < n on S A to 
write dVL\ = d6 2 + sin 2 9dQ 2 and €4 = sin 3 9d9 A £3 with e n being the volume form of the 
unit S n . Our D5-brane is wrapping S 3 at a constant 9 whose value will be determined 
dynamically by a non-zero k Fl flux along x 12 on its world-volume. We choose the gauge 
such that C RR is smooth around 9 = to have 



c rr = (2*ls) 3 N / e _ 1 cog3 ^ _ 2 \ ^ A £3 (3 3Q) 

V4 V 3 3 / 

With these data, it is rather straightforward to compute the above D5 action, except a 
subtle point that the k Fl charge along x 12 is represented by an imaginary F 12 = iF [5]. 
One way of seeing this is to consider a Wick-rotated Lorenzian situation with x\ being 
the time direction, where the Fl string charge would be unambiguously given by a real 
F\2- Upon going back to our Euclidean situation by the Wick-rotation of xi, the Fl flux 
will transform to an imaginary value. 

The resulting D5-brane action density on x% 2 after integrating over S 3 and £4 is 



or ' - N sin 3 9 



- F> + - (cosO - -cosH - -) F , (3.31) 



where (U, 9) is yet undetermined position of the _D5-brane, and 

4 5 vr 6 iV 2 T 6 
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The k Fl charge on its world-volume is now identified as a conserved flux 

/ „ \ 



k 



5sf 
5F 



N 

T 



F sin 3 9 



3 cos 9 + cos 3 9 + 2 



F 2 



(3.33) 



/ 



which can be solved for F in terms of (U,9), and the effective Hamiltonian density one 
obtains after a Legendre transform with the conserved flux F becomes 



h 



,™ + kF='- 



N 



u - 3 
C <TT T 



\ 



sin 6 + (3 cos0 - cos 3 9 - 2 + ^ 



(3.34) 



and we have to minimize this with respect to (U, 9) to find the fc-wall tension in our 
framework. It is trivial to see that the D5-brane sits at the IR end U = Ut, while the size 
of S 3 inside S A given by the polar angle 9 must be determined by solving the following 
equation 



sin 



(3.35) 



9_\ _ k_ 
2)~N 

As the function of 9 in the left-hand side is a monotonic function in the range [0, ir] with 
values between (0, 1), one checks that k (N — k) corresponds to 9 <-» (tt — 9) in the 
solution. With the solution of 9 and U = Ut, the fc-wall tension is finally given by 

- 3 T 3 32tt 3 T 3 



N 



7T 



T^^v^sin 2 ^ o _ 



N 2 sin 2 9- 



M, 



KK 



27 



-k(N-k) 



M, 



KK 



(3.36) 



with the desired property of = T^-k- Note that we obtain the precise Casimir scaling 
Tk ~ k(N — k), contrary to the N = 4 SYM result in Ref. [5]. For a small 4 1, one 



also confirms that the result reduces to the tension of k D2-branes (13.251) . 



32vr 3 , Ar T 3 
-kN- 



27 



M 



KK 



k i 



(3.37) 



4 Fate of Zn vacua with fundamental flavors 



In this section, we come to our second motivation, that is, studying what happens to the 
Zj<j- vacua with a presence of fundamental flavors in a probe approximation, or equivalently 
in a quenched approximation, via gauge/gravity correspondence. We first consider the 
system of D3/D7-branes that is dual to iV = 4 SYM perturbed by a small number of 
N = 2 fundamental hypermultiplets represented by Nf probe D7-branes in a iV D3- 
brane background [6]. We are especially interested in the dependence on the mass of the 
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flavors given by the asymptotic distance between D7 and D3-branes. We also calculate 
the energy lift of Z^-vacua in the massless case analytically, and finally perform a similar 
analysis in the model by Sakai-Sugimoto [7J for a QCD-like theory with massless chiral 
quarks. 



Although the appearance of discrete Zn vacua in the deconefined phase geometry ( 12. 8ft 
of TV = 4 SYM is a sub- leading effect of large iV limit, one can still advocate a quenched 
approximation where one neglects back-reaction of the probe D7-branes to the N = 
4 SYM dynamics including the mechanism of selecting Zn vacua. We point out that 
this is just one type of approximation which is not solely based on the large N limit, 
because back-reactions of D7-branes may well be of the same sub-leading order of the 
previous mechanism of Z^ vacua in gravity. However, we feel that the present quenched 
approximation is a useful one to pursue due to its practical calculability and the possibility 
of its comparison to lattice QCD in the same quenched approximation. 

Working in the quenched approximation in the gravity side means that we take the 
previous Z N vacua of 



as a given background and consider the probe dynamics of D7-branes embedded in it. 
As we neglect possible back-reactions, questions regarding the fate of Z N vacua become 
those of the probe _D7-branes, such as whether the energy of /c'th vacuum is lifted or not. 
We can address this question by studying the energy of the probe D7-brane in the fc'th 
gravity vacuum (14.381) . 

The action of a probe D7-brane is 



where the Chern-Simons term is irrelevant for our purpose. One should study the equa- 
tions of motion of the embedding X M (£) and the world- volume gauge field F = dA in the 
presence of the background B NS in ( 14.381) . and in general these two are coupled to each 
other, except the case of trivial vacuum k = where one can turn off F = consistently. 
This special case was analyzed in Ref.[23J with varying hypermultiplet mass m q , where 
they found a phase transition of meson-melting at T c ~ m q . Our problem is a more 
sophisticated version of theirs. 



4.1 D3/D7 system 




(4.38) 




(4.39) 
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To parameterize D7-brane embedding in the background of (12.81) . it is convenient to 
introduce the rectangular coordinate x 4 ' 5 ' 6,7 ' 8,9 such that r and S 5 coordinates are related 
to them by 



dr 2 + r 2 dn 2 5 = J2(d. 



IX 



i\2 



(4.40) 



i=4 
?6 



and in fact they are nothing but the original R coordinates transverse to the D3-branes. 
Asymptotically at r —* oo, our D7-brane spans four flat directions out of x 4-9 , say x 4-7 , 
being a point at x 8 ' 9 with a distance 2^l 2 m q from the center. This fixes the UV boundary 
condition specified by the mass of the fundamental flavor m q . As r goes to the infrared 
region, the D7-brane would feel an attraction toward the horizon at r = Tn and it bends. 
One easily finds that the spherical symmetry on x 4-7 is a symmetry of the situation, and 



the D7-brane wraps S -sphere in x given a point in x ' . Moreover D7-brane would 
move only along one axis on the x 8 ' 9 -plane, and one can take x 9 = without loss of 
generality. Then, the embedding is simply given by a map between x 8 and the radius p 

(x 8 ) 2 + (x 9 ) 2 + p 2 and 



of S 3 inside x 4 7 . One notes that r 2 



dr 2 + r 2 dtt\ 



(dx 8 ) 2 + (dx 9 ) 2 + dp 2 + p 2 dQ 2 3 



(4.41) 



so that one can equally describe the embedding by a map between x 8 and r, which will 
be chosen from now on. We will choose the world-volume coordinates £ a by (ig, Xi, r, f^), 
and x 8 is a function of r specifying the embedding shape that is determined dynamically 
with the specified UV boundary condition at r — > oo, 



x 8 (r — ► oo) = 2nl 2 m q 



(4.42) 



as we discuss in the above. The induced metric g* ab on the D7 world- volume is then 



d? 2 

ab D7 



+ 



r 

Z 2 

L 2 



7r 4 T 4 L S 



dt p — )— o?x,' 



/ 8/ \dx s (r) 

[r - x 8 (r)-^ 
r 2 — (x 8 (r)) 2 




+ 



|x 8 (r) 

7T 4 T 4 L 8 
.4 _ ^4^8 



dQ 2 3 



dr 2 



(4.43) 



from which one easily computes the action density along Xi,2,3 after integrating over 
and S 3 , 



-,D7 



dr 



(x 8 (r)) 



. 3 
2\ 2 



A+(B» s + 2nP s F 0r y 



(4.44) 
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with 

Wi-ggn \tzf^f + (m)\ , (445) 



r 4 



r-* 8 (r)^) ^(r)V , tt 4 T 4 L 8 



r 2 _ ( x 8( r ))2 \ dr j r 4 - vr 4 T 4 L 8 



where we turn on F 0r along (£e, r) as a possible world- volume flux in response to the given 
background -Ev s with 



2tt/2 , rH 



We are assuming symmetry along the thermal circle without loss of generality. 

The 
motion, 



The equation of motion of the gauge field is simple to solve; ^f- is a constant of 



— ^- ~ . ' v = r^C , (4.47) 

SFor y/A + (^ + 2nP s F 0r f 

with some constant C and we factored out r\ for later convenience. One can solve the 
above for {B^ r s + 2irl 2 s F 0r ^j, and the remaining dynamics of the D7 embedding x 8 (r) will 
then be described by the Routhian, 

R ° = ^~I dr (4 < s + ^) S = sfej / * We -<*•«>')' , 

(4.48) 

given a value of C. The quest is to scan possible range of C and the resulting £>7-brane 
embedding shape from solving the above R^ 7 , to minimize the original action density 
s^ 7 . Once the final configuration is found, one can compute its energy density from s® 7 
to address the question of energy lift of vacua. The necessary numerical analysis, as 
one varies the quark mass m q , is beyond the scope of the present work. Identifying a 
phase transition would be an interesting future direction. 

However, we make several important remarks. Note that in our coordinate parametriza- 
tion of D7 embedding, there is an inequality 

x 8 {r) < r , (4.49) 

and the equality x 8 (r) = r happens precisely when S 3 shrinks to zero, where the D7-brane 
stops going into a lower r, that is, it is the end point of the D7 embedding in the radial 
direction r. If this happens at r = r > r# in the solution, the D7-brane doesn't touch 
the horizon, and we are in the phase of un-melted mesons in the gauge theory side. It is 
important to observe that in this phase the D7-brane covers only part of the r) cigar 
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D of the background (12. 8p . that is, a cylinder topology of [r , oo] Xitg inside the cigar D. 
Looking back the expression (I4.47j) . one readily finds that the constant of motion C must 
be zero in this phase considering the point r = r , and this in turn says that 

+ = , (4.50) 

in all region of the D7-brane world-volume. Then the analysis of the embedding shape 
and the resulting energy density from s^ 7 becomes precisely equal to the case of trivial 
vacuum k = or without B N /F fluxes at all. One concludes that Z^- vacua are not 
lifted and persist to exist in this phase in quenched approximation. One interprets that 
the world-volume flux F nullifies the background B NS flux to minimize s® 7 . In the gauge 
theory point of view, as the world-volume gauge symmetry on the D7-brane is the flavor 
global symmetry in the gauge theory side, this points out an interesting phenomenon of 
a dynamically generated monodromy of flavor symmetry along the thermal circle that 
counteracts the existing Z^ Wilson line in the Z^-vacua, to make the Zjv-vacua viable 
even in the presence of fundamental flavors. One mathematical remark is that (I4.50P is 
possible precisely because the topology of D7-brane in r) is cylindrical, so that there 
is no topological restriction of its J F value. 

In the other phase where x 8 (r) < r remains true until it hits the horizon r = r#, the 
Z)7-brane meets the horizon and the mesons are melted. A crucial point in this phase is 
that the D7-brane now wraps the whole cigar D of (t_E,r) with topology of a two-plane, 
and there is an important topological restriction 

/ F — 2nm , me Z . (4.51) 
Jd 

Suppose one tries to put C = 0, or equivalently ^^B^. 3 + F 0r = 0, to minimize the action 
s® 7 as before, then 

/ F = -^[ B NS = - 2 -^f2,Z , (4.52) 

except the trivial vacuum k = 0. In other words, in the /c'th vacuum it is not possible to 
completely nullify the Zn Wilson line B N due to a topological restriction, and one has a 
non-vanishing C ^ or ^If -^0r S + Fq t ^ in the non-trivial Z N vacua. For each integer 
m, C must be determined to satisfy (14.51 j) . and one should choose m that minimizes the 
energy s^ 7 . Looking at the expression (14.441) of s^ 7 , one sees that the resulting energy is 
always higher than the trivial vacuum k = where ^j-^or 5 + For = 0, and we conclude 
that the Z^ vacua are lifted in the phase of melted mesons. 
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The question at which phase the system finds its solution given the UV boundary 
condition of the flavor mass m q should be addressed numerically, but for a sufficiently 
large ^> 1 it is natural to expect the former phase where the mesons are un-melted and 
the vacua survive, while in the extreme case of m q = it is obvious that x 8 (r) = and 
we are in the latter phase of melted mesons with Z^-vacua lifted. One naturally expect 
a phase transition to happen at some point T c ~ m q . An interesting question is whether 
the critical temperature T c depends on k or not because of the additional contribution 
from the B N /F fluxes to s^ 7 , whose answer is beyond the present paper. 

In the simplest case of m q = where one can consistently set x 8 (r) = to be in the 
melted-meson phase, it is possible to compute the energy lift of the Z^-vacua analytically. 
Note that we have A = 1 in (I4.44p with x 8 (r) = 0, and the equation of motion for the 
gauge field (I4.47j) simplifies as 



r 3 (B» s + 2rrl 2 s F 0r 



(B» s + 2nl 2 s F 0r ) 



= r%C , (4.53) 



which can be solved for F, 



Or 



2tt/ : 



; B 0r + . ° = . (4.54) 



As discussed before, C can not take any continuous values, but is restricted to only 
discrete values to satisfy the topological constraint (14.51 p . Explicitly, integrating both 
sides of (14.541) over D, one has 

2nk C/3 f°° , 1 2nk Cf3r H /■«> 1 

27rrra = + / dr . = + — / dr - , (4.55) 

N 2nP s Jr H Uj_y_ C 2 N 2nP s Ji Vr 6 - C 2 ' 1 J 

so that given (m, k), the C is determined by 

C Tdr J „ =^2F 1 (llic 2 )=2ir(m+±-)^L= , (4.56) 



Vf 6 - C 2 2 V3'2'3' J \ N J V^IiV 

where we have used Th = 7rTL 2 = ttT (AnggN) 1 1 2 . Once C is determined as a function 



of (m, k), it is straightforward to calculate the energy density s|? 



s ■ 

D7 



Q D7 m 



(3r% r°° , f f 6 A irg s N 2 T 3 r°° IA / f 6 , 
dr = — r = ■ / dr = — r 



k =o 2 6 ir 5 g s l 8 s Ji \^f 6 - C 2 J 4 A \^f 6 - C 2 

(4.57) 
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where we have subtracted the value of the trivial vacuum k = to see the energy lift of 
the (m, k) state. For a fixed k, one should choose m that minimizes the energy to finally 
find the energy lift of the fc-vacuum, 

A e , = min mex (s D E % m k) -s°\ =o ) . (4.58) 

From f )4.57p . one has to minimize C 2 , and an inspection of ( 14. 56ft shows that this is 
achieved by the smallest value of \m+ 4|. Therefore, our result is invariant under k — > 
k + N and k —>■ (N — k) , which are the desired properties due to Zjq nature and charge 
conjugation. For low lying fc-vacua with ^ < 1, we have m = to minimize the energy, 
and an easy calculation gives us 

Ae^vr^T 3 , ±<1 . (4.59) 

One also finds that Ae^ has a discontinuous slope at jq = \- 

4.2 The model of Sakai and Sugimoto 

The Sakai-Sugimoto model is obtained by introducing probe D8/D8 branes in the Wit- 
ten's ZM-brane background, that span S 4 x R 3,1 x U and are point-like in a^-direction. In 
the original weak-coupling DA/ D8/ D8-brane picture, one finds left(right)-handed chiral 
quarks from D4 — D8(DA — D8) string spectrum, and at low energy the resulting 4D 
gauge theory flows precisely to QCD with massless chiral flavor quarks. One therefore ex- 
pects the gravity picture to capture some of the interesting dynamics of this semi-realistic 
gauge theory, at least for low enough energy regime. In the confined geometry (I3.15P where 
(x4,U > Ukk) makes a cigar shape with vanishing a; 4 at U = Ukk, the D8/D8-bra.ne 
pair must join with each other at IR, realizing spontaneous chiral symmetry breaking in 
an intuitively geometric way. Our current interest is however on the opposite case of the 
deconfined geometry (I3.19p . where remains finite everywhere, and we instead have a 
cigar shape D with (t^, U > Ut)- In this case, each D8 and Z)8-brane separately meets 
the horizon at U = Ut and their world-volume embedding is simply parameterized by a 
constant £4 position. We will consider a single pair of D8/D8-branes for simplicity. 

As we consider D8/D8-branes in the Z N - vacua of the background, 

1 f B NS = ^rdUB^ = 2 -^ , (4.60) 



2-kIIJd 2rrZf Ju T uu N 

one can easily find that the embedding shape can still be consistently given by a constant 
£ 4 from the equation of motion, and we will take this without further detail. One can 
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also be convinced by translational symmetry along X4. Because the Chern-Simons term 
of the D8/D8 doesn't play any role here, the actions of D8 and D8 are identical, and the 
total action is two times the DBI action of the D8-brane, 

S E = 2/i 8 J d 9 £ e^det (g* + B NS * + 2nP s F) , (4.61) 

where we choose the world-volume coordinates £ a to be (£e, Xi, U, It is straightforward 
to compute the above, turning on a possible world- volume gauge field Fqu along the cigar 
D of (tE, U). After integrating over S 4 x ts, the action density in the spatial indirection 
becomes 

1.7-1 

N 2 r°° b / o 

s e = 5 i g / dtf ^ V 1 + W + 2vr/ s 2 F 0l ;) 2 . (4.62) 

3-2 A n-giL 2 T JUt 

The equation of motion for F ;y is simply -j^p- = const., which we parameterize by 



8F 0U 

Ui (B™ + 2irP s F 0U _ 
'ou) 2 



= U^C , (4.63) 



y/l + (B$f + 2nP s F , 
with a dimensionless constant C, and from this one easily solves for Fqu as 

1 „ 1 C 



ou 



(*) - 



As in the previous case of D3/D7 system, an important point is a topological restric- 
tion 

p /*oo 

F = /3 dUF 0U = 2nm , m e Z , (4.65) 

JD JUt 

and this determines the constant C in terms of given (k, m) as follows, 



C dU ; = — 2 Fi — , -, — , C 2 = m + — — -i = m + — 



^^5 _ C 2 3 HlO^'lO' y V /3C/ T V N)Ang s NlsT 



2C „ / 3 1 13 „ Q \ / fc \ 4tt 2 / 2 / fc 

J 4^ 

(4.66) 

where 2 -Pi is the hypergeometric function. Once C is given, one can easily calculate the 
difference in the action density se from the k — vacuum where (^B^f + 2ti1 2 s Fqij^ = 0, 
and after sorting out the coefficients, the final answer is 



se - s E 




[ dUU% \ ^— 1 . (4.67) 



For a given /c-vacuum, one should find the integer m such that it minimizes the above 
result, and that is the final energy lift of the /c'th vacuum Ae^. Because C 2 is a monotonic 
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increasing function of 



m + 



A' 



and this again gives us Ae 
i 

2 



one should take the smallest possible value of m + 
Aek+N and Ae^ = Ae7v_ fe . The slope is discontinuous in 



N 



k — ^tfc+V 

27 1 



the middle jr = \. For small ^ < 1, one has m = and one can easily find that 



C 



k 



8vr g s Nl s T N 



and 



Ae k « Q^g s l s k 2 T A = 2 3 3 4 vr 2 (^ M iV 



k 

N« 1 



k 2 T 4 



(4.68) 



(4.69) 



where we have used <7 S Z S 
system (14.591) . 



2wM KK 



-. This is qualitatively different from the case of D3/D7 



5 Conclusion 

In this work, we first study fc-walls with k ~ in the Witten's D4-brane background of 
pure SU(N) Yang-Mills theory in large N limit, largely motivated by a recent work in 
AdSz, x S 5 by Armoni, Kumar and Ridgway [5]. We propose and check consistency of 
the picture that k (Euclidean) D2-branes blow up into a A^S^-brane wrapping 5" 3 inside 
S A of the background via Myers effect. We compute the tension of A;- wall by performing 
suitable U-duality to transform the A^SS-brane into a D5-brane. Our result is a precise 
Casimir scaling behavior of the /c-wall tension ~ k(N — k). We note that the same 
Casimir scaling was also obtained for the fc-quark flux tube tension in Ref . [24] . 

Our second subject is to consider fate of Z N -vacua in the presence of fundamental 
flavors in quenched approximation via gauge/gravity correspondence. In our study of 
D3/D7 system, we point out an interesting phenomenon of phase transition as one varies 
^r, which separates the phase of un-melted mesons with Z/y- vacua survived and the phase 
of melted mesons with ZAr-vacua lifted. An interesting question is whether the critical 
temperature depends on k, which should be addressed by a numerical analysis in the 
future. 

In the special case of m q = 0, we calculate the energy lift of fc'th vacua analytically. 
Our result is consistent with the ^Ar-nature and the charge conjugation where k — > (N—k). 
We also find a tower of stable states parameterized by integers m above the lowest energy 
state. We do the similar analysis in the Sakai-Sugimoto model for an approximate QCD- 
like theory. 
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